Statistical transmutation in doped quantum dimer models. 
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We prove a "statistical transmutation" symmetry of doped quantum dimer models on the square, 
triangular and kagome lattices: the energy spectrum is invariant under a simultaneous change of 
statistics (i.e. bosonic into fermionic or vice-versa) of the holes and of the signs of all the dimer 
resonance loops. This exact transformation enables to define duality equivalence between doped 
quantum dimer Hamiltonians, and provides the analytic framework to analyze dynamical statistical 
transmutations. We investigate numerically the doping of the triangular quantum dimer model, 
with special focus on the topological Z2 dimer liquid. Doping leads to four (instead of two for 
the square lattice) inequivalent families of Hamiltonians. Competition between phase separation, 
superfluidity, supersolidity and fermionic phases is investigated in the four families. 



The discovery of exotic liquids such as topological Z 2 
spin liquids (SL) [10 is one of the challenges of modern 
condensed matter physics. Anderson proposed that the 
parent (insulating) state of the high-temperature super- 
conductors is in fact a SL, the resonating valence bond 
(RVB) state and that spin-charge separation (and su- 
perconductivity) will occur under doping (j^: the origi- 
nal electron fractionalizes into two emergent particles, a 
holon carrying the charge quantum and a spinon carry- 
ing the spin quantum. Although the original electron is 
a fcrmion, there has been a long-standing debate regard- 
ing the actual statistics of holons and spinons in such 
a "deconfinement" scenario. In this context, Rokhsar 
and Kivelson Q introduced the quantum dimer model 
(QDM) as a simple effective model to describe magnet- 
ically disordered phases. The basis assumption here is 
that the spinon spectrum is gapped [H[ (strictly speak- 
ing, in the QDM it is infinite) and the dimers between 
nearest neighbor (NN) sites mimic fluctuating SU(2) sin- 
glets of paired electrons. The underlying microscopic ex- 
change interaction leads to effective attraction/repulsion 
between dimers and dimer-flips along closed loops (see 
below). Doped QDM, where dimers (i.e. pairs of elec- 
trons) are removed from the system, leading to itiner- 
ant holes, have also been studied 0, 041] ■ There, holes 
and dimers are strongly coupled due to hard-core con- 
straints. Variants of doped QDM have been also con- 
structed to physically describe polarized spinons induced 
by a magnetic field @. Naively, one expects holons to 
be of fermionic nature (while spinons should be bosonic) 
but, in spite of this naive expectation, the statistics of 
holons has been debated over the last 20 years. Earlier 
work suggested that holon excitations were fermions flp| . 
It was then argued that, in fact, holon statistics is dic- 
tated by energetic considerations and under some condi- 
tions a holon could become a bosonic composite through 
binding of a flux quantum [8(, the "vison" in the QDM 
language. This was indeed observed recently by exact di- 
agonalization techniques in the square lattice Q where, 
by varying the ratio of dimer kinetic energy vs holon 




FIG. 1: (color online) Dimer coverings corresponding to dif- 
ferent lattices. Light blue areas represent flippable plaquettes 
and arrows represent the allowed single hole hopping process, 
i) Square lattice: The hole hopping is between sites in the 
same sublattice. ii) Triangular lattice: A hole can hop in the 
different directions schematized by the arrows, iii) kagome 
lattice: The light blue area shows one of the possible flip- 
pable plaquettes of length 10. 



kinetic energy, one can transit from a regime where low- 
energy quasi-particles behave as fermions to a regime 
where they behave as bosons. 

In this letter we prove an exact duality transforma- 
tion that provides a framework to analyze this dynam- 
ical statistical transmutation. More precisely, we study 
the generic Hamiltonians of QDM's doped with holes de- 
picted in Fig. [U which contain three contributions, Hj 
and Hy respectively the flipping and potential energy 
terms for the dimers of amplitudes J and V and H t the 
hopping term for vacancies (named "holes" ) of amplitude 
t. The Hilbert space corresponds to all dimer coverings 
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with a fixed hole density. The potential term is propor- 
tional to the number of flippablc plaqucttcs (which is not 
conserved by dimer flips and hole motion). Special atten- 
tion is devoted to frustrated lattices, the triangular lat- 
tice with edge-sharing triangular units and the kagome 
lattice with corner- sharing triangular units, which we 
compare to the results obtained for the square lattice Q ■ 
The exact symmetry of these Hamiltonians states that 
their spectrum is invariant under simultaneous transfor- 
mations of the statistics of the holes, bosons into fermions 
or vice-versa, B J 7 , and of the sign of (all) dimer ki- 
netic amplitudes, J <H> — J, V and t being unchanged. 

Equivalence classes - For fixed values of the mag- 
nitudes of </, V and t, one can define eight different 
Hamiltonians by changing the signs of J and/or t and 
by choosing bosonic or fermionic statistics for the mo- 
bile holes. Using the exact symmetry mentioned above 
and proven below, these eight different Hamiltonians can 
be grouped in pairs with identical spectra, hence defin- 
ing four non-equivalent families. This is summarized in 
Fig. [2] More precisely, the established equivalence is 
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FIG. 2: (color online) Equivalence classes ("models") and re- 
lations between the eight Hamiltonians introduced in the text. 
Starting form the unfrustrated (J > 0, t > 0) bosonic Hi 
Hamiltonian, one can combine any change of statistics and/or 
change of signs of J and/or sign of t to obtain H2, ■ ■ ■ , H$. 
Equivalence relations are shown by arrows: full dark grey ar- 
rows (dotted light grey arrows) are valid for all lattices (for 
vacancies hopping on the square lattice only). 



the square lattice by H 1 -f-» Hi O H3 O H4, on one 
hand, and by H§ -H- Hq O H7 <R- H$ on the other hand, 
were dubbed as "Perron-Frobenius" and "non-Perron- 
Frobenius" Hamiltonians and studied in Q while, for the 
triangular lattice, only the unfrustrated Hi Hamiltonian 
(representing family (a)) was studied in Q. 

As the transformation we shall establish below change 
bosons into fermions (and vice versa), we can start by 
assuming, without lost of generality, that the (bare) 
holons in the doped QDM are bosons. We implement a 
two dimensional Jordan- Wigner transformation on these 
bosons to change their statistics to fermionic. In con- 
trast to one dimensional systems, this resulting trans- 
formed Hamiltonian is highly non-local [0 - 14 1 requiring, 
in general, a mean field approximation to proceed fur- 
ther, at least from an analytical point of view. Next, we 
show that the non-local terms can be absorbed by using 
a different representation for the dimer operators, which 
keeps their bosonic character. This key feature provides 
then an elegant proof of the "statistical transmutation 
symmetry" of these models. 

Proof of statistical transmutation symmetry - It is 
convenient to write the Hamiltonian in a second quan- 
tized form, by introducing creation operators b\ ■ for a 
dimer siting between sites i and j and holes by oper- 
ators at. In our conventions, dimers between sites i 
and j are created by spatially symmetric operators b\ ■, 
both operators (b\j and a\) are bosonic and mutually 
commuting. It is instructive to notice that the dimer 
bosonic operators can be thought as bilinears of "elec- 



trons" operators: b 
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). In terms 



of these operators we implement the hard-core constraint 
a l a i+J2z "i i+zbi,i+z = lj where the sum runs over NN of 
site i. Let us call V the projector on the subspace where 
the constraint holds. In the following we use systemat- 
ically the projected Hamiltonian VHV (= HV since all 
terms of H commute with V), which prove to be very 
useful later. 

Let us now apply 2-D Jordan- Wigner (JW) transfor- 
mation on the holon operators (l2j : 



a, = e"** 



fi 



(1) 



between Hi and Hi, both represented by the family la- 
bel (a) and in the same way for the pairs of Hamiltoni- 
ans i?2i-i and Hn with i = 2,3,4 corresponding respec- 
tively to families (b) , (c) and (d) . We argue that such an 
equivalence is a general feature of doped-two-dimensional 
QDMs, independently of the lattice geometry. Physically, 
classes (a) and (b) ((c) and (d)) describe systems with 
doped polarized spinons (doped holons), and are called 
"spinonic" QDMs ("electronic" QDMs). Note that, for 
the square (or hexagonal) lattice, an additional t -O- — t 
symmetry reduces the number of non-equivalent fami- 
lies from 4 to 2. In Fig. [3J the two families defined for 



where <fo = fj fj ar si T j ~ n) and 73 = Xj + iyj 

is the complex coordinate of the j-th hole. Using that 
arg(Ti — Tj) = arg(rj — r») ± tt it follows immediately 
that two /-operators in different sites anticommute and 
the constraint and the phase <pi can be written equally 
in terms of fi or a, operators. 

In order to understand the consequences of transforma- 
tion |T]) on the Hamiltonian, the hopping of holons can 
be written, for an arbitrary lattice, as a sum of three-site 
Hamiltonians 
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Making use of the transformation (QJ we obtain 
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In the last equation we have changed boson operators a,i 
by fermionic ones fi at the cost of introducing non-local 
interactions. In other words we have written a boson as 
a composite particle consisting of an electron with an at- 
tached flux. We can use the same procedure to change the 
dimer operators from the fermionic to the bosonic repre- 
sentation: in the fermionic representation of the dimers, 
we attach the remaining flux to the fermions to obtain 
a bosonic representation. In order to define new dimer 
operators including the phases 4> in their definition, we 
must be able to write <f> in terms of operators bij. This 
can be performed in the following way: When applied 
into the subspace projected by V, we can change <f> by </> 
in the exponentials, where 
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The remaining non-local exponential operators are re- 
written in terms of <p which can be absorbed by defining 
bi.j = e 1 ^" 1 "^' bij. It is a simple matter to see that 
operators b\j are also bosonic, as it should be to make 
sense as dimer operators. Then, transformation (jXJ) to- 
gether with the definition of bij allow us to change the 
statistics of holes from bosonic to fermionic. After this 
transformation the hopping term is written in terms of 
operators bij and fi and the hopping amplitude changes 
to t = ie^+argfo-TO-argCT,--!-*)^ Let us now investi- 
gate the effect of this transformation specifically for each 
lattice : 

i) Square lattice. The Hamitonian is defined by 
Hi = Hj + Hy + H t , where: 
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Before going to the details of the Jordan- Wigner trans- 
formed Hamiltonian let us mention some well known 
gauge transformations which can be performed. For zero 
doping, a simple gauge transformation on the dimers can 
be done to show the equivalence of the Hamiltonians with 
J and — J. In the case of non-zero doping, changing the 
sign of the holes wave function with a wave vector trans- 
formed as k — > k + (-7T, 7r) one obtains easily the equiva- 
lence between Hamiltonians with t and — t. We can now 



proceed to the computation of the Jordan- Wigner trans- 
formed Hamiltonian. The transformed Hamiltonian be- 
comes, after some algebra, H 2 = Hj^-j + Hy + H t , 
where the tilde means that dimers and holes are created 
by operators b^ and f* respectively. In the dimer kinetic 
Hamiltonian the amplitude J has been changed by — J. 

ii) Triangular lattice. - Here the Hamiltonian Hi is 
written using: 
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and similar expressions for and Ha for the other ori- 
entations of the rhombi [Hj]. corresponds to the 
holon hopping on the up-triangles, the total hopping 
Hamiltonian being H t = H^ + H^\ It is important 
to note that, on such non-bipartite lattices, the equiva- 
lence t <-> — t does not hold anymore, as verified numer- 
ically. This implies that there are now four families of 
equivalent Hamiltonians, instead of two for the square 
case. On the triangular lattice, the dual Hamiltonian 



reads Ho 



H 



(J)^(-J) 



H K P + with the same 

convention as before. A corollary of this result is that, as 
for the square lattice, the Hamiltonians with J and — J 
are equivalent in the zero doping case. 

in) Kagome lattice. - On this geometry, all dimer 
resonant loops (of length a = 6, 8, 10, 12) involving a 
single hexagon become important and should be consid- 
ered [3]. Nearest-neighbor hole hopping as depicted in 
Fi g, rr jean be introduced as in [l7| (the reader can refer 
to 16|, [l7[ for more details) . The JW transformed Hamil- 



tonian is obtained in the same way as before and since 
it is interesting on its own, it will be presented elsewhere 
[l8| . Here, we only summarize the main results. By us- 
ing the JW transformation followed by (more involved) 
gauge transformations on dimer and hole operators, [19( 
one can write the dual Hamiltonian with all kinetic dimer 
amplitudes J a changed into — J a . As before, dimers and 
holons are created by ¥ and p respectively. 

Statistical transmutation and choice of Hamiltonian - 
The above statistical transmutation symmetry is of great 
importance in analytic and numerical investigations of 
doped QDMs, establishing their phase diagrams. In 0], it 
was shown that, in the square lattice, by varying the ratio 
J 1 1 or the doping, the system undergoes a series of phase 
transitions. One of those phase transitions corresponds 
to a dynamical transmutation, between a phase where el- 
ementary low-energy quasi-particles are fermionic holes, 
to a phase where the low-energy quasi-particles acquire 
a bosonic nature. It was in fact predicted that holon and 
vison could pair up, leading to such statistical transmu- 
tation Q . Then, if the microscopic Hamiltonain is chosen 
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FIG. 3: (color online) Qualitative phase diagrams of the 
four families of doped QDM's (a-c) as labeled in Fig. [2] on 
the triangular lattice versus doping (x) and Vy|J| for fixed 
t/J — 0.5. Inset: Energy vs normalized magnetic flux (£) 
corresponding to family (c) for fixed V/\J\ = 0.5, t/J — 0.5 
and x — 0.25. The E vs £ plots corresponding to families (a), 
(b) and (d) are qualitatively similar. 



such that holons are fermions, the fermionic phase corre- 
sponds to a weak-coupling regime and the bosonic phase 
to a strong-coupling one. This picture gets interchanged 
if holons in the microscopic Hamiltonian are chosen to 
be bosons. Hence, using this duality equivalence one 
can always choose the most relevant microscopic QDM 
Hamiltonian to be in a weak-coupling regime (depending 
on the point in the phase diagram). 

Phase diagrams - We now complement the exact re- 
sults with a numerical study of the phase diagram of the 
four non-equivalent families of Hamiltonians in the tri- 
angular lattice where a topological (Z2) liquid can be 
stabilized at zero doping |2j, |20j. The statistics of the 
dressed excitations is studied using the method devel- 
oped in Q by investigating the node content of the wave 
functions. Family (a) in Fig. [3] is the only unfrustratcd 
case and was studied with Green Function Monte Carlo 
methods in Q. In this model bare holons arc bosonic 
and remain representative of the physical excitations in 
the entire region of the phase diagram, as happen also 
in the Perron-Frobenius square lattice version, studied 
in 0. The unfrustrated model shows a superfluid phase 
in a wide range of the phase diagram. As the effective 
charge is Q = 2e even in the vicinity of the Z2 dimcr 
liquid region, the composite particle is not produced by 
a simple confinement as in the VBS region 0, [2l[ (see 
below). The situation is even more interesting when J 
is changed into — J, or cquivalcntly, bosons are changed 
into fermions (families (c) and (d) in Fig. [3]) . Most of the 



phase diagram is occupied by what we dubbed a "com- 
plex phase" as the statistics of dressed excitations does 
not correspond solely to bosons or fermions. Similarly to 
the square lattice, we find a phase separated region at 
low doping. Family (d) also shows, in addition, an inter- 
esting fermion reconstruction at large doping, which we 
call "fermionic phase" . 

In order to provide a better understanding of the con- 
ducting phases of the triangular lattice, we insert in 
the torus an Aharonov-Bohm flux of strength cf> = £0o 
with < £ < 1 and where (f>o = hc/e the elemen- 
tary magnetic flux. As done in Ref. @ on the square 
lattice, one can reduce the finite-size effects by consid- 
ering arbitrary boundary conditions in the y-direction. 
All this is achieved by the Peierls substitution t —> 
cxp(2iTr((;u x + 2ri/\/3uy) ■ Ui/L)t with Ui the hopping 
direction on the triangular lattice of linear size L. A 
superfluid would be characterized by well-defined min- 
ima in the ground state energy separated by a finite bar- 
rier in the thermodynamic limit. A contrario, a typical 
signature of (weakly interacting) fermions, a fiat energy 
profile is expected even on such a small cluster 22j . Here, 
we report that the ground state energy has well-defined 
minima quantized at half a flux quantum for all family 
of models at x ~ 0.25, compatible with a charge Q = 2e 
superconductor (see inset of Fig. 3(c)). 

Discussion and perspectives - In this paper we have 
shown that the nature of the excitations in doped QDM 
is a much more subtle question than what one would 
naively think. We have rigorously established equiva- 
lence classes between QDM Hamiltonians (see Fig. [2]) 
which provides a powerful tool to identify the nature of 
the dressed hole excitations. We claim that this dual- 
ity relation is a generic feature of QDM independently 
of the details of the lattice. From a dynamical point of 
view, we provide evidence of statistical transmutation on 
a frustrated lattice, analogously to the case of the square 
lattice. In some families of Hamiltonians we have found 
a complex phase in which dressed excitations can not be 
understood in terms of solely fermionic or solely bosonic 
degrees of freedom. Nevertheless, all these models reveal 
flux quantization in units of half a flux quantum, consis- 
tent with the idea that the dimer background leads to 
effective particles of charge 2e. 
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